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(2)  Objectives: 


The  objectives  of  the  present  research  include  development  of  a  suitable  CFD  solver  for  solving  incompressible  flow 
past  airfoils  suitable  for  flying  wing  micro  air  vehicle  (MAV)  design.  Suitable  airfoils  would  have  to  be  selected  for 
the  purpose.  This  would  be  followed  by  development  of  a  structural  analysis  solver  that  would  structurally  model 
the  airfoil  and  integrate  this  solver  with  the  flow  solver  in  order  to  perform  the  fluid  structure  interaction  studies  for 
flexible  airfoils.  Compute  flow  past  selected  flexible  airfoils  and  compare  their  aerodynamic  and  structural 
performance  with  corresponding  rigid  airfoils.  Implementation  of  smart  concept  in  the  structure  solver  in  order  to 
improve  the  flexible  airfoil  performance.  Perform  the  fluid  structure  interaction  studies  by  incorporating  smart  active 
layer  and  explore  whether  such  an  embedded  layer  can  improve  the  aerodynamic  and/or  structural  behaviour  of  the 
flexible  airfoil.  The  studies  are  to  be  performed  over  a  moderate  range  of  Reynolds  number  and  a  fairly  wide  range 
of  angles  of  attack  to  cover  a  reasonably  wide  variety  of  flow  conditions.  The  scope  of  the  project  also  includes  flow 
solution  past  three  dimensional  finite  wings  with  some  optimization  of  the  wing  planform. 


(3)  Status  of  effort: 

To  achieve  the  objectives  of  the  present  project,  the  following  milestones  were  achieved  sequentially. 

•  An  existing  two  dimensional  structured  grid  based  finite  volume  Navier  Stokes  flow  solver  was  modified 
into  an  unstructured  grid  based  flow  solver.  It  is  capable  of  solving  flow  past  complicated  multi  body 
geometries  over  low  to  moderate  Reynolds  numbers.  It  was  validated  on  a  few  benchmark  problems 
successfully  before  being  applied  for  the  present  research  work. 

•  The  flow  solver  was  used  to  perform  flow  simulations  on  a  large  number  of  reflexed  airfoils,  mainly  Eppler 
series  airfoils,  which  are  candidate  airfoils  for  flying  wing  MAV  design.  The  simulations  were  performed 
over  moderate  range  of  Reynolds  number  and  angle  of  attack.  The  airfoils  were  assumed  to  be  rigid. 

•  A  finite  element  based  Structural  Analysis  Solver  was  developed  which  would  model  the  structural 
behavior  of  the  airfoil.  The  solver  was  validated  on  several  benchmark  problems. 

•  The  smart  layer  comprising  of  embedded  piezoelectric  sensors  and  actuators  was  incorporated  in  the 
structure  solver  and  was  validated  on  several  benchmark  problems. 

•  The  fluid  and  structure  solvers  were  interfaced  to  compute  the  fluid  structure  interaction  problem  which 
occurs  when  considering  an  elastic  wing.  Implementation  of  smart  concept  in  the  coupled  fluid- structure 
solver  to  enhance  the  aerodynamic  and  structural  characteristics  of  the  airfoil  depending  on  the  properties 
of  the  embedded  piezoelectric  sensors  and  actuators. 

•  Carrying  out  fairly  large  number  of  flow  simulations  by  considering  both  rigid  and  flexible  airfoils  under 
uniform  and  gusty  flow  conditions.  Impact  of  a  sinusoidal  gust  flow  was  studied  for  two  different  airfoils 
because  MAVs  are  likely  to  encounter  gusts  during  flight  and  their  aerodynamic  and  structural  performance 
may  be  severely  affected  in  such  unsteady  environment.  For  gusty  flow  simulations  two  airfoils  were 
considered,  namely,  Eppler  340  (E340)  and  a  3%  uniform  thickness  airfoil  with  same  mean  camberline  as 
E340. 

Certain  observations: 

From  the  flow  simulation  results  that  were  undertaken  in  the  present  project  it  was  obvious  that  the  airfoil  shape 
would  need  to  be  optimized  to  obtain  optimum  aerodynamic  performance  over  a  wide  range  of  angle  of  attack, 
instead  of  using  existing  airfoil  shapes.  Flowever,  such  optimization  would  demand  lot  more  effort  and  time  which 
could  affect  achievement  of  other  broad  goals  of  the  present  research.  Therefore  it  was  identified  as  a  study  which 
could  be  undertaken  in  future.  The  existing  two  dimensional  structured  incompressible  Navier  Stokes  solver  was 
extended  to  two  dimensional  unstructured  solver  in  order  to  attain  desired  mesh  fineness  around  airfoil  with  overall 
lower  number  of  cells  in  the  computational  domain  as  compared  to  structured  grid  case  and  also  to  attain  the 
flexibility  of  using  more  complex  body  geometries,  including  multi  element  airfoils.  Flowever,  extension  of  this 
solver  to  three  dimensional  unstructured  NS  solver  could  not  be  successfully  achieved  because  of  considerable 
complexities  which  were  encountered  and  which  could  not  be  resolved  within  the  time  frame  of  the  project.  This 
solver  would  be  developed  over  the  next  several  months  so  that  it  could  be  used  in  future  to  solve  flow  past  finite 
MAV  wing.  Investigators  would  like  to  work  on  this  development  for  some  more  time  so  that  fluid  structure 
interaction  studies  may  later  be  extended  to  full  three  dimensional  flexible  wing  geometries.  Also  an  attempt  could 
be  made  in  optimizing  location  of  smart  patches  over  such  a  three  dimensional  wing  planform. 


(4)  Abstract: 


A  flying  wing  MAV  configuration  has  several  advantages  like  reduced  weight,  cleaner  aerodynamics,  lesser  number 
of  control  surfaces  and  actuators,  easier  attainment  of  stringent  size  restrictions  as  compared  to  a  conventional 
configuration  with  tail.  However,  a  tailless  wing  may  not  be  stable.  Therefore,  the  pitching  moment  of  the  airfoil 
should  be  low  and  preferably  positive.  This  could  be  achieved  by  reflexed  airfoils  like  Eppler  and  S  series  airfoils  or 
low  pitching  moment  airfoils  like  MH  series  airfoils.  A  thin  wing  could  also  be  considered  which  has  the  same  mean 
camber  line  as  that  of  a  conventional  airfoil,  including  the  trailing  region  reflex  for  the  same  purpose.  Thin  wings 
could  be  more  suitable  to  facilitate  active  control  and  thereby  improve  the  aerodynamic  and  structural  characteristics 
of  the  airfoil  for  widely  varying  flow  conditions.  The  airfoil  could  be  actively  controlled  by  suitable  sensors  and 
actuators.  The  present  paper  focuses  on  numerical  investigation  of  the  aerodynamic  and  fluid  structure  interaction 
aspects  of  a  few  reflexed  and  non-reflexed  Eppler  airfoils  as  well  as  thin  airfoil  of  constant  3%  thickness  and 
identical  mean  camberline  as  that  of  a  reflexed  airfoil,  namely  E340  airfoil.  These  airfoils  could  all  be  potential 
candidates  for  a  flying  wing  MAV  application.  The  structural  simulation  performed  by  the  finite  element  method  is 
coupled  with  a  finite  volume  flow  solver  via  a  partitioned  solution  approach.  Most  published  approaches  use  a 
partitioned  solution  scheme,  where  specialized  programs  for  fluid  and  structural  dynamics  are  loosely  coupled.  This 
approach  basically  transfers  in  each  time  step  surface  pressure  and  shear  force  from  the  fluid  to  the  surface  of  the 
structure  and  updates  the  geometry  of  the  fluid  domain  according  to  the  computed  structural  displacement.  A 
piezoelectric  sensor-actuator  layer  is  embedded  within  the  airfoil  in  order  to  improve  its  aerodynamic  and  structural 
response  over  a  wide  range  of  angle  of  attack.  Both  a  uniform  flow  at  various  angles  of  incidence  as  well  as  cosine 
wave  gust  has  been  used  to  study  the  steady  as  well  as  unsteady  characteristics  of  the  airfoil  over  a  modest  range  of 
low  to  moderately  high  Reynolds  number  ( 1 0  to  8X104)  which  is  compatible  with  practical  MAV  applications.  The 
fluid  and  structure  solvers  have  been  validated  independently  on  certain  benchmark  problems.  Subsequently  they 
have  been  integrated  together  for  performing  the  fluid  structure  interaction  studies.  For  performing  the  flexible 
airfoil  studies,  the  three  different  flexible  cases,  namely,  (i)  undamped  uncontrolled  case  (ii)  damped  case,  with 
natural  damping  of  2%,  and  (iii)  active  control  case,  i.e.,  with  AFC  actuator  and  PVDF  sensor  have  been  considered. 
The  aerodynamic  characteristics  of  rigid  airfoils  as  well  as  flexible  airfoils  have  been  compared.  The  active  control 
seems  to  be  most  impressive  in  terms  of  its  vibration  damping  capabilities. 
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(10)  Archival  Documentation: 

A  MAV  flies  in  a  low  Reynolds  number  regime  due  to  its  low  flight  speed  and  small  dimensions.  Such  flow  is  often 
accompanied  by  laminar  boundary  layer  separation,  formation  of  laminar  separation  bubble,  transition  and 
accompanying  low  lift-to-drag  ratio.  The  fluctuations  in  wind  speed,  which  can  be  comparable  to  MAV  flight  speed, 
make  both  the  instantaneous  flight  Reynolds  number  and  angle  of  attack  vary  substantially  and  results  in  an 
unsteady  flight  environment.  In  the  present  work  a  few  Eppler  airfoils  have  been  chosen  for  investigation.  They 
include  Eppler  61  (E61),  Eppler  330  (E330),  Eppler  334  (E334)  and  Eppler  340  (E340)  respectively.  E61  is  a  non- 
reflexed  airfoil  whereas  the  remaining  are  reflexed  airfoils.  Reflex  angle  (5)  is  the  negative  camber  angle  found  at 
the  trailing  edge  of  reflexed  airfoils.  The  reflex  angle  of  E330,  E334  and  E340  are  5=5.3°,  2.6°  and  6.6°  respectively. 
Uniform  flow  past  these  airfoils  have  been  computed  for  Reynolds  number  (Re)=103,  104,  3X104,  4.6X104  and  8x  104 
respectively  at  various  angles  of  attack.  The  essential  features  of  the  flowfield  have  been  captured.  Pressure  contours, 
streamlines,  time  averaged  airfoil  surface  pressure  distribution,  variation  of  pitching  moment  with  angle  of  attack 
and  drag  polar  plots  have  been  included.  Uniform  flow  as  well  as  gusty  flow  has  been  solved  for  E340  airfoil  at 
Re=8x104  at  various  angles  of  attack  by  treating  the  airfoil  initially  as  rigid  and  later  as  flexible.  Similar  flow 
calculations  have  been  performed  for  a  3%  thick  airfoil  having  same  mean  camber  line  as  E340  airfoil  (henceforth 
referred  to  as  3%  constant  thickness  airfoil)  at  Re=8x104  at  various  angles  of  attack.  The  flexible  airfoil  is 
constrained  at  30%  chord  for  structural  calculations.  The  flow  field  characteristics  and  time  history  of  force 
coefficients  for  both  rigid  as  well  as  flexible  airfoil  cases  have  been  investigated.  Time  history  of  the  displacement 
of  the  trailing  edge  of  the  two  airfoils  were  obtained  (Figure.  15  &16).  Fast  Fourier  Transform  of  the  trailing  edge 
displacement  is  used  to  determine  the  frequency  of  vibration  of  the  airfoil.  In  the  present  work  a  loosely  coupled 
approach  has  been  followed  for  carrying  out  the  fluid  structure  interaction  studies  using  two  different  fluid  and 
structure  solvers  which  exchange  information  back  and  forth. 


The  Numerical  Scheme 
The  Flow  Solver 


In  order  to  solve  the  flow  field  two  dimensional  Navier-Stokes  equations  for  incompressible  flow  in  Cartesian 
coordinates  are  used,  which  can  be  written  using  indicial  notation  as  follows: 


<9u 


J-  =  0 


(1) 

<3P  ,  dxn  (2) 

The  above  governing  equations  are  solved  in  their  non  dimensional  form  in  the  present  flow  solver.  The  constitutive 


<5u  d  /  \  dp  dij. 

-  +  -  U:U:  =-  —  +  - ~ 

dt  dxj v  ‘  ’  dx;  dXj 


relation  between  stress  and  strain  rate  for  a  Newtonian  fluid  is  used  to  link  the  components  of  the  stress  tensor  T-  to 

velocity  gradients.  Earlier,  a  finite  volume  flow  solver  was  developed  based  on  ‘Consistent  Flux  Reconstruction’ 
(CFR)  scheme  of  Roy  and  Bandyopadhyay  [1]  which  was  originally  developed  on  a  structured  grid  (referred  earlier 
as  existing  flow  solver).  Pressure  Poisson  equation  is  used  for  computing  the  pressure  field  which  gives  a  divergence 
free  velocity  field.  In  the  present  work  an  unstructured  triangular  mesh  based  ‘CFRUNS’  solver  has  been  developed 
by  Harichandan  and  Roy  [2]  which  is  based  on  the  original  CFR  scheme.  This  solver  is  used  to  compute  the  flow 
field  around  the  airfoil.  Details  of  the  solver  are  available  in  the  above  reference.  The  initial  unstructured  mesh  is 
generated  using  GAMBIT8.  Subsequently  when  the  flexible  airfoil  shape  changes,  the  flow  domain  is  re-meshed 
using  a  suitable  moving  grid  or  grid  re-meshing  procedure  based  on  spring  analogy  or  Master/Slave  concept  (Lian  et 
al  [3]).  Since  the  airfoil  is  flexible  it  would  deform  under  the  action  of  air  pressure  and  shear  stress  acting  on  its 
surface  from  time  to  time.  Since  the  solver  is  based  on  moving  grid  approach,  geometric  conservation  law  is 


enforced  to  ensure  consistent  flow  calculations.  A  suitable  upwinding  scheme  is  incorporated  in  the  solver  for 
enhancing  its  stability  at  higher  Reynolds  numbers  (Wu  and  Hu  [4]).  The  size  of  the  computational  domain  in  terms 
of  airfoil  chord  length  (c)  are  as  follows:  length=40c,  with  10c  distance  upstream  and  30c  distance  downstream  of 
airfoil;  height=20c,  which  is  equally  disposed  on  top  and  bottom  portions  of  the  airfoil.  Approximately  35000  cells 
have  been  used  to  discretise  the  computational  domain.  Grid  convergence  studies  have  been  performed  for  the 
CFRUNS  solver  which  are  reported  in  Harichandan  and  Roy  [2].  The  following  inlet  velocity  distribution  has  been 
used  to  create  gusty  flow  past  E340  airfoil. 

U=U0 

V=V0  y>4,  y<-4  (3) 

and  V=  Vo  +  a  cos(cot)  -4  <  y  <  4 

where  U0,  Vo,  a  and  co  are  1.0,  0,  0.8  and  ji/4  respectively  and  y  is  the  distance  measured  in  the  direction  normal  to 
the  flow  direction  with  y=0  located  at  airfoil  mid  chord.  The  time  period  of  a  gust  cycle  is  8  non  dimensional  time 
units. 

The  Structure  Solver 

Smart  materials  could  be  surface  bonded  or  embedded  as  sensor/actuator  layers  in  the  host  structure  to  monitor/ 
control  the  behavior  of  the  structure.  An  Active  Fiber  Composite  (AFC)  smart  layer  which  is  embedded  within  an 
airfoil  could  be  used  to  monitor  and  actively  control  its  shape  and  thereby  help  in  enhancing  its  aerodynamic 
performance  under  varying  operating  conditions.  AFC  consists  of  unidirectional  PZT  fibres  embedded  into  epoxy 
matrix  and  a  separate  electrode  layer  with  an  inter-digital  pattern  (Bent  et  al.  [5];  Bent  and  Hagood  [6]).  In  the 
present  investigation  in  order  to  solve  for  the  structural  response  of  the  smart  airfoil  to  the  applied  aerodynamic 
loads  a  finite  element  (FE)  model  has  been  developed.  An  eight  noded  two  dimensional  quadratic  isoparametric 
element  with  5  mechanical  degrees  of  freedom  and  1  electrical  degree  of  freedom  per  node  is  used  to  formulate  the 
FE  model.  The  behaviour  of  smart  structure  is  governed  by  the  equations 

A  numerical  code  has  been  developed  based  on  the  formulated  FE  model.  The  variable  thickness  of  the  airfoil  is 
modeled  as  a  beam  of  variable  thickness.  According  to  the  thickness  of  the  airfoil  at  any  chordwise  position  the 
element  stiffness  is  calculated  separately  and  finally  the  global  stiffness  matrix  is  assembled.  In  the  present  case  a 
graphite/epoxy  laminated  airfoil  is  considered.  At  each  calculation  step  the  FE  code  provides  the  displacement  data 
for  the  airfoil  surface  nodes  to  the  fluid  code.  The  FE  code  can  run  both  for  uncontrolled  as  well  as  actively 
controlled  conditions.  The  actively  controlled  airfoil  results  have  been  obtained  based  on  a  constant  voltage 
application.  Three  different  flexible  cases  have  been  considered  in  the  present  investigation,  namely,  (i)  undamped 
uncontrolled  case  (ii)  damped  case,  with  natural  damping  of  2%,  and  (iii)  active  control  case,  i.e.,  with  AFC 
actuator  and  PVDF  sensor.  Based  on  the  displacement  data  the  fluid  domain  is  re-meshed.  A  moving  grid  approach 
is  used  for  this  purpose.  The  fluid  code  is  run  till  steady  state  results  are  obtained  for  the  new  shape  of  the  airfoil. 
The  steady  state  body  surface  pressure  and  shear  stress  distribution  is  fed  back  to  the  FE  code  to  obtain  fresh 
displacement  data.  This  interactive  fluid-structure  calculations  are  performed  till  suitable  convergence  criteria  are 
satisfied.  The  time  stepping  for  the  fluid  and  structure  solvers  have  been  judiciously  chosen  to  suit  the  characteristic 
timescale  of  relevant  physical  phenomena. 

Results  and  Discussion 

Results  of  the  investigations  concerning  both  rigid  and  flexible  airfoil  cases  are  discussed  below.  Figure  1  shows  the 
(a)  instantaneous  pressure  contours  and  (b)  streamlines  of  flow  past  E340  airfoil  at  angle  of  attack  (a)=12°, 
(Re)=103.  The  flow  separates  from  the  leading  edge  of  the  airfoil  and  a  Von  Karman  vortex  street  is  formed  in  the 
wake  of  the  airfoil.  At  this  Reynolds  number  the  flow  is  expected  to  undergo  a  laminar  separation  from  the  leading 
edge.  The  Reynolds  number  is  not  high  enough  to  trigger  a  reattachment  which  will  lead  to  formation  of  a  laminar 


(4) 

(5) 


separation  bubble  near  the  leading  edge.  Due  to  negative  camber  near  the  trailing  edge  of  the  airfoil,  the  cross-over 
of  pressure  distribution  occurs  on  the  airfoil  surface  at  moderate  angles  of  attack.  This  is  clearly  visible  from  Figure 
2  at  a=8°.  This  crossover  phenomenon  creates  negative  lift  over  portion  of  the  airfoil  near  trailing  edge  thereby 
reducing  the  overall  airfoil  lift.  However,  it  is  favorable  from  pitching  moment  point  of  view.  A  fairly  extensive 
study  has  been  carried  out  on  reflexed  Eppler  airfoils  (E325-E339)  as  a  part  of  the  present  research  which  reveals  the 
need  for  choosing  the  proper  reflex  angle  to  get  the  correct  choice  of  maximum  lift  and  pitching  moment  (Roy  [7]). 
Increased  reflex  decreases  the  maximum  lift  coefficient  and  increases  the  pitching  moment  and  vice  versa.  Figure 
3(a)  and  (b)  show  the  pitching  moment  about  quarter  chord  and  drag  polar  respectively  of  E330  and  E334  airfoils 
which  clearly  indicates  the  effect  of  variation  of  reflex  angle.  Similar  studies  were  also  performed  on  a  number  of 
reflexed  S  series  airfoils  which  indicated  comparable  findings.  E61  is  a  non  reflexed  airfoil  as  mentioned  before. 
The  instantaneous  streamlines  past  E61  are  shown  in  Figure  4  for  Re=4.6x  1 04at  a=T2°,  15°  and  25°  respectively.  A 
laminar  separation  bubble  is  observed  near  the  trailing  edge  at  a  =12°  which  gradually  seems  to  enlarge  with 
increase  in  angle  of  attack.  At  a=25°,  the  bubble  has  ultimately  led  way  to  a  large  scale  separation  and  consequent 
vortex  shedding.  The  time  averaged  pressure  distribution  on  the  airfoil  for  various  angles  of  attack  are  shown  in 
Figure  5.  The  results  compare  reasonably  with  those  reported  in  Savaliya  et  al  [8].  All  the  above  studies  were 
performed  assuming  a  rigid  airfoil. 

For  the  FE  solver,  simulations  have  been  initially  carried  out  for  2-D  simply  supported  rectangular  plates  and  the 
results  have  been  validated  for  static  and  dynamic  conditions  with  those  available  from  literature.  Using  the  above 
code,  behavior  of  cross  ply  and  angle  ply  laminates  has  been  investigated.  Different  loadings  like  uniformly 
distributed  load,  impulse  load  and  transverse  harmonic  loads  have  been  tested.  Uncontrolled  and  active  controlled 
conditions  have  been  simulated.  The  active  control  achieved  with  embedded  AFC  layer,  as  expected,  gives  much 
more  favorable  characteristics  than  the  uncontrolled  case.  Figure  6  is  a  schematic  representation  of  a  rectangular 
composite  plate  with  AFC  actuator  layer  and  PVDF  sensor  layer.  The  response  of  the  plate  under  uncontrolled  and 
active  control  for  impulse  load  and  transverse  harmonic  load  is  shown  in  Figures  7(a)  and  (b)  respectively  (Kerur 
and  Ghosh  [9]).  Figure  8  shows  the  airfoil  laminate  arrangement  with  AFC  actuator  &  PVDF  sensor. 

Figure  9  shows  the  instantaneous  streamlines  and  vorticity  contours  which  show  the  alternate  formation  of 
separation  bubble  on  upper  and  lower  surfaces  of  flexible  E340  airfoil  under  gust  for  Re=8x104,  a=0°.  The  flow 
history  is  represented  from  t=32  to  40  time  units  at  interval  of  one  time  unit.  Careful  inspection  of  the  flow  field 
pictures  at  t=32  and  40  reveals  the  fact  that  gust  cycle  is  completed  in  8  time  units.  Figure  10  shows  instantaneous 
streamlines  for  gust  flow  past  flexible  3%  thick  airfoil  at  non  dimensional  times  t=24.0  to  32.0,  at  interval  of  one 
time  unit,  at  Re=3x104,  a=12°.  A  laminar  separation  bubble  visible  on  the  rear  bottom  side  of  the  airfoil  at  t=24.0 
vanishes  at  t=25.0  due  to  increase  in  angle  of  attack.  A  large  separation  bubble  develops  over  the  upper  surface  of 
the  airfoil  at  t=26.0  which  enlarges  at  t=27.0  due  to  further  increase  of  angle  of  attack  but  is  still  attached  to  the 
airfoil  surface.  At  t=28.0  a  large  eddy  structure  leaves  the  airfoil  surface  and  forms  secondary  eddy  structures  near 
leading  and  trailing  edges.  At  t=29.0  the  flow  over  the  upper  surface  reattaches  due  to  decrease  of  angle  of  attack 
and  the  shed  vortex  gets  convected  into  the  wake.  At  t=30.0  an  almost  streamlined  flow  is  observed  over  the  airfoil. 
At  t=31.0  a  laminar  separation  bubble  forms  near  the  leading  edge  on  the  bottom  surface  of  the  airfoil.  This  bubble 
moves  to  the  rear  portion  of  the  lower  surface  at  t=32.0.  Figure  11  shows  lift  and  drag  coefficient  history  for  E340 
airfoil  in  gust  flow  at  Re=8x104,  a=0°.  It  is  noticed  that  since  the  inlet  is  located  10  chord  lengths  upstream  of  the 
leading  edge  of  the  airfoil  it  takes  10  non  dimensional  time  units  for  it  to  impinge  on  the  airfoil.  As  a  consequence, 
the  variations  in  lift  and  drag  coefficients  start  growing  at  around  that  time.  The  lift  force  has  a  time  averaged  value 
of  zero  as  a=0°.  Figure  12  shows  the  lift  and  drag  coefficient  history  for  3%  airfoil  in  gust  flow  at  Re=3x104,  a=12°. 
Figure  13  shows  a  comparative  representation  of  time  varying  lift  acting  on  3%  airfoil  in  gust  flow  at  Re=3x104, 
a=12°  for  the  flexible  and  rigid  cases.  It  is  obvious  from  the  figure  that  the  three  different  flexible  cases  show  very 
small  difference  between  each  other  in  terms  of  their  effect  on  lifting  characteristic  of  the  airfoil.  However,  it  is 
noticed  from  the  figure  that  there  is  some  difference  with  the  corresponding  rigid  airfoil  case.  This  makes  it  obvious 
that  the  aerodynamic  characteristics  of  the  airfoil  does  get  affected  due  to  its  flexibility.  The  extent  of  variation 
would  depend  on  the  material  properties  of  the  flexible  material  and  the  consequent  stiffness  of  the  section.  Figure  14 
shows  a  zoomed  view  of  the  lift  variation  produced  on  the  E340  airfoil  for  the  three  different  flexible  cases.  Figure 
15  and  16  show  the  non  dimensional  deflection  of  the  trailing  edge  of  the  airfoil  for  two  different  cases.  The 
displacement  of  the  trailing  edge  is  expressed  as  percentage  of  the  maximum  thickness  of  the  airfoil.  It  is  found  that 
the  active  case  is  most  effective  in  dampening  in  the  vibrations,  more  radically  for  the  3%  airfoil  case.  Figures  17 
and  18  show  the  frequency  of  vibration  of  the  trailing  edge  of  the  E340  and  3%  thick  airfoils  respectively.  The 
frequency  spectra  are  obtained  by  using  Fast  Fourier  Transform  of  the  displacement  time  history  of  training  edge 


node.  The  comparatively  higher  frequency  of  the  E340  airfoil  could  be  linked  with  its  higher  stiffness  due  to  higher 
average  section  thickness. 


Conclusions 

A  two  dimensional  unstructured  grid  based  incompressible  Navier  Stokes  solver  ‘CFRUNS’  was  developed  to 
perform  the  present  computations  for  flow  past  various  reflexed  and  non-reflexed  Eppler  airfoils,  namely,  E61,  E330, 
E334  and  E340.  The  initial  unstructured  grid  to  discretise  the  fluid  domain  was  generated  using  GAMBIT®  software. 
The  fluid  domain  was  remeshed  whenever  necessary  by  using  a  moving  grid  technique  accompanied  by 
implementation  of  geometric  conservation  law  in  the  solver.  Validation  of  the  ‘CFRUNS’  solver  has  been  reported 
in  Flarichandan  and  Roy  [2].  The  airfoil  flow  computations  were  performed  over  a  range  of  Reynolds  numbers  103, 
104,  4.6xl04  and  8X104  for  various  angles  of  attack.  In  the  above  simulations  the  airfoils  were  assumed  to  be  rigid. 
Both  effect  of  uniform  flow  as  well  as  gusty  flow  past  various  airfoils  was  investigated.  Instantaneous  pressure 
contours,  streamlines,  time  averaged  surface  pressure  distribution,  force  and  moment  coefficients  etc  were  used  for 
performing  the  necessary  analysis.  The  essential  features  of  the  flowfield  were  captured  through  the  flow 
simulations.  A  finite  element  based  smart  structure  solver  was  developed  as  a  part  of  the  present  research.  The 
solver  was  validated  by  applying  it  on  some  benchmark  problems  which  are  reported  in  the  paper.  A  loosely  coupled 
fluid  structure  solver  was  subsequently  developed  and  applied  to  solve  the  fluid  structure  interaction  problem  arising 
in  flexible  airfoils.  Uniform  flow  as  well  as  gusty  flow  was  solved  for  flexible  E340  airfoil  at  Re=8x104  for  various 
angles  of  attack  as  well  as  3%  thick  airfoil  at  Re=3x104.  Lift  time  history  have  been  compared  between  the  rigid  and 
flexible  cases.  All  the  flexible  cases  produce  similar  lift  characteristics  but  differ  to  some  extent  with  the  rigid  airfoil 
values.  The  airfoil  trailing  edge  displacement  time  history  reveals  that  the  time  varying  loading  on  the  airfoil  in  a 
gusty  situation  produces  oscillations  of  the  airfoil  which  could  be  most  effectively  dampened  by  the  embedded 
active  piezo  layer.  The  differing  trailing  edge  vibration  frequencies  could  be  associated  with  the  varying  stiffness  of 
the  E340  and  3%  thick  airfoils.  Finally,  it  was  evident  from  the  present  study  that  development  of  a  successful  fluid 
structure  interaction  model  initiates  the  scope  of  broader  study  related  to  vibration  suppression  using  tailored  AFC 
layers  to  the  shape  control  of  airfoil  to  improve  the  in-flight  performance  of  a  MAY. 
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Figure  1:  Instantaneous  (a)  pressure  contours  and  (b)  streamlines  for  flow  past  E340  airfoil  at  a=12°,  Re=103 


(C) 

Figure  4:  Instantaneous  streamlines  for  uniform  flow  past  rigid  E61  airfoil  at  Re=4.6*104,  a=(a)  12°  (b)  15° 

and  (c)  25°  respectively 


Figure  5:  Time  averaged  pressure  distribution  on  E61  airfoil  for  uniform  flow  at  Re=4.6*104  for  various 

angles  of  attack 


Figure  6:  Plate  Configuration  with  AFC  Actuator  and  PVDF  Sensor 
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Figure  7:  Uncontrolled  and  active  controlled  non-dimensionaiised  centre  deflection  of 
(PVDF/0/90/0/90/AFC)  cross  ply  laminate  subjected  to  (a)  impulse  load  (b)  transverse  harmonic  load 


Figure  8.  Airfoil  laminate  arrangement  with  AFC  actuator  &  PVDF  sensor 


_t  =  33 


Figure9:Instantaneous  streamlines  and  pressure  contours  for  gust  flow  past  flexible  E340  airfoil  at  non 
dimensional  time  t=32.0  to  40.0,  at  intervals  of  one  time  unit,  at  Re=8xl04,  a=0° 


t=  27 


t=  29 


Figure  10:  Instantaneous  streamlines  for  gust  flow  past  flexible  3%  thick  airfoil  having  same  mean  camber 
line  as  E340  airfoil  at  non  dimensional  times  t=24.0  to  32.0,  at  interval  of  one  time  unit,  at  Re=3xl04,  a=12° 
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Figure  11  Lift  and  drag  coefficient  history  for  E340  airfoil  in  gust  flow  at  Re=8xl04,  a=0° 
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Figure  12  Lift  and  drag  coefficient  history  for  3%  airfoil  in  gust  flow  at  Re=3x104,  «=1 2' 


Non-dimensional  Time 


Figure  13  Time  varying  lift  acting  on  3%  airfoil  in  gust  flow  at  Re=3xl04,  a=12°  for  the  flexible  and  rigid 

cases. 
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Figure  14  Lift  variation  produced  on  the  E340  airfoil  for  the  three  different  flexible  cases 
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Figure  15  Trailing  edge  tip  deflection  for  Eppler  340  airfoil  for  the  three  flexible  cases 
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Figure  16  Trailing  edge  tip  deflection  for  Eppler  340  airfoil  for  the  three  flexible  cases 


Figure  17  Vibration  frequency  spectra  for  Eppler  340  airfoil  for  the  three  flexible  cases,  Re=8xl04,  a=0° 


Figure  18  Vibration  frequency  spectra  for  3%  thick  airfoil  for  the  three  flexible  cases,  Re=3xl04,  a=0' 


